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INTRODUCTION 


A number  of  azide  compounds  including  lead  azide  are  known  to 
detonate  when  subjected  to  voltages  or  fields  of  sufficient  strength 
(Ref  1,  4).  A considerable  reduction  in  the  values  of  voltage  or 
field  necessary  for  initiation  upon  application  of  low- intensity  ir- 
radiation of  appropriate  wavelength  was  recently  reported  by  us 
(Ref  5).  The  mechanisms  involved  in  these  phenomena  are  not  under- 
stood in  detail  although  they  appear  to  involve  electronic  processes 
(Ref  5).  A complete  experimental  characterization  of  these  phenomena 
is  essential  to  their  detailed  understanding.  We  have  completed  the 
first  step  in  such  a characterization  by  studying  the  response  of 
unilluminated  lead  azide  single  crystals  to  electric  fields  using 
evaporated  Au- contacts  in  a particular  geometric  configuration. 

The  results  are  of  interest  in  another  important  context.  The 
problem  of  the  sensitivity  of  lead  azide  to  static  electricity  and 
to  electric  discharge  has  received  considerable  attention  because  of 
practical  considerations  relating  to  safety  of  storage  and  handling 
(Ref  6,  8).  Studies  in  this  area  have  typically  been  made  on  lead 
azide  powders.  It  would  appear,  however,  that  a fundamental  under- 
standing of  this  sensitivity  must  be  based  upon  a knowledge  of  the 
effects  of  electric  fields  on  lead  azide  in  its  pure  crystalline 
form. 


The  scope  of  our  study  was  limited  somewhat  by  the  fact  that  our 
experiments  involve  the  destruction  of  samples,  and  lead  azide  crystals 
suitable  for  our  purposes  are  not  routinely  available  in  quantity. 
Hence,  we  concentrated  on  two  goals.  The  first  was  to  determine 
whether  or  not  reproducible  threshold  values  of  either  voltage  or 
field  are  responsible  for  initiation.  The  second  was  to  study  the 
current- vo It age  characteristics  of  voltages  below  and  near  the  initia- 
tion condition.  Deviations  from  linearity  can  reveal  electronic  ef- 
fects such  as  charge  injection  which  may  be  involved  in  causing  the 
initiation  of  lead  azide. 


EXPERIMENTAL 


Single  crystals  of  a- lead  azide  were  grown  by  slow  cooling  from 
an  ammonium  acetate  solution,  following  the  method  developed  by 
Garrett  (Ref  9)  in  this  laboratory. 
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Several  crystals  were  wafered  using  a string  saw  with  a weak 
solution  of  ceric- ammonium-nitrate  wetting  the  string.  These  wafers 
were  then  polished  on  both  sides  using  very  fine  aluminum  oxide 
abrasives.  It  should  be  noted  that  the  wafers  were  not  all  sliced 
perpendicular  to  the  same  crystallographic  direction. 

Vacuum-evaporated  gold  electrodes  were  applied  to  both  sides  of 
the  wafers.  These  gold  electrodes  were  then  contacted  by  gold  wires 
using  dots  of  silver  paint.  Figure  1 shows  a typical  wafer  with  the 
gold  electrode  applied.  Sample  thicknesses  varied  from  .019  to 
.076  cm  in  this  series  of  experiments.  The  electrode  area  was  kept 
constant  at  .04  cm2.  The  crystal  mounting  is  shown  schematically  in 
Figure  2. 

A sandwich  type  electrode  configuration  was  chosen  because  sur- 
face currents  can  then  be  essentially  eliminated;  the  entire  sample- 
electrode  system  was  potted  in  RTV,  a polymeric  electrode  potting 
compound  which  further  aids  in  prohibiting  the  flow  of  current  by  any 
path  other  than  through  the  bulk  of  the  sample.  By  keeping  the 
evaporated  electrode  diameter  large  compared  to  the  sample  thickness, 
the  results  can  be  analyzed  on  the  basis  of  one-dimensional  planar 
current  flow. 

The  circuitry  and  experimental  procedures  were  kept  as  simple  as 
possible.  The  circuitry,  shown  in  Figure  3,  consists  of  a regulated 
D.C.  power  supply,  picoammeter,  and  strip  chart  recorder.  RG-59/U 
shielded  cable  was  used  to  connect  the  components.  A 2-ma  fuse  was 
placed  in  series  with  the  sample  to  protect  the  input  resistor  of  the 
picoammeter  from  the  large  transient  current  which  occurs  at  deto- 
nation. 

The  sample  chamber,  shown  in  Figure  4,  was  evacuated  by  a fore- 
pump with  an  in-line  liquid  nitrogen  trap.  With  samples  of  the  size 
used  in  these  experiments,  the  chamber  is  not  damaged  by  the  explo- 
sion. Our  experience  has  been  that  performing  the  experiments  under 
this  type  of  vacuum  is  sufficient  to  eliminate  any  dependence  of  the 
results  on  atmospheric  humidity  conditions. 

A serious  procedural  difficulty  arose  because  of  the  possibility 
that  lead  azide  crystals  subject  to  high  electric  fields  can  experi- 
ence localized  decomposition  or  other  damage.  This  would  act  as  a 
memory  of  prior  treatment  and  might  be  expected  to  influence  both 
conductivity  and  critical  field  for  detonation.  One  way  around  this 
difficulty  would  involve  using  a fresh  crystal  of  particular  thick- 
ness for  each  separate  voltage  measurement.  This  procedure  would, 
however,  have  demanded  far  more  samples  of  accurately  controlled 
thickness  than  we  could  prepare. 
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To  minimize  possible  deviations  caused  by  memory  effects,  we 
decided  to  subject  each  sample  to  the  same  voltage  history.  Our  pro- 
cedure involved  one -minute- on,  one-minute-off  application  of  voltage, 
increasing  the  voltage  in  100V  increments  as  illustrated  in  Figure  5. 
The  procedure  was  continued  until  the  sample  detonated  with  current 
being  monitored  continuously.  Because  sample  thickness  was  a variable, 
it  should  be  noted  that  all  samples  were  not  subjected  to  the  same 
field  history. 

The  experiments  were  performed  at  room  temperature  (,v300K) , and 
the  samples  were  not  illuminated. 


RESULTS 


Figure  5 shows  the  typical  current  response  to  the  applied  volt- 
age steps.  There  is  an  undershoot  when  the  voltage  is  turned  off 
which  is  not  shown.  Seven  samples  of  varying  thickness  were  investi- 
gated in  this  manner.  The  magnitude  of  the  current  I just  prior  to 
detonation  and  the  average  electric  field  strength  at  detonation 
Vdet/k  are  s^own  Table  1.  It  should  be  noted  that  the  detonation 

does  not  necessarily  occur  immediately  upon  application  of  the  volt- 
age. The  time  lag  was  observed  to  vary.  We  have  also  plotted  the 
logarithm  of  the  current  observed  after  one  minute,  versus  the  loga- 
rithm of  the  applied  voltage.  These  results  appear  in  Figures  6-12. 


DISCUSSION 


Table  1 shows  that  the  voltages  at  which  detonation  occurred 
(Vdet)  ranged  from  800  to  2700V.  Thus  it  would  hardly  make  sense  to 

speak  of  a threshold  voltage.  One  would  expect  a voltage  threshold 
if  the  voltage  drop  occurred  primarily  at  the  electrodes.  If  the 
voltage  drop  across  the  sample  is  relatively  uniform,  the  bulk  of  the 
sample  experiences  an  electric  field  near  the  average  field,  defined 
as  V^et  divided  by  sample  thickness  L.  The  average  fields  at  which 

detonation  occurred  are  rather  more  closely  grouped,  these  ranging 
from  30,250  to  42,000  V/cm.  It  thus  appears  that  there  is  a thres- 
hold average  field  at  which  detonation  occurs  in  this  field  range. 

The  scatter  in  values  would  then  be  accounted  for  by  sample  differ- 
ences beyond  our  control,  including  perhaps  deviations  in  the  uni- 
formity of  voltage  drop  or  differences  due  to  crystal  orientation. 
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Note  that  the  threshold  field  values  for  the  four  thinnest  samples 
(.019  - .025  cm)  are  all  higher  than  those  for  the  thicker  samples 
(.043  - .076  cm).  This  would  be  expected  if  damage  resulting  from 
exposure  to  fields  tends  to  lower  the  field  threshold,  since  the 
thicker  samples  experienced  longer  cumulative  exposure  times  to 
fields. 

We  now  turn  to  a discussion  of  the  current-voltage  characteris- 
tics. Since  no  attempt  was  made  to  eliminate  capacitance  from  the 
external  circuit,  the  initial  capacitive  spike  which  occurs  imme- 
diately upon  application  of  the  voltage  does  not  accurately  represent 
current  through  the  sample.  (The  peak  current  was  linear  with  volt- 
age over  the  entire  voltage  range  for  all  samples.)  The  current  at 
the  end  of  one  minute,  although  not  strictly  steady  state,  was  taken 
to  represent  the  response  of  the  current  through  the  crystal  to 
applied  voltage. 

The  log  I- log  V plots  of  the  seven  samples  are  not  all  similar 
(Figures  6-12).  Portions  of  some  of  them  show  regions  where  current 
depends  rather  clearly  on  V,  V2,  and/or  V3  and  some  contain  discon- 
tinuities. The  variance  in  these  I-V  characteristics  is  rather  puzz- 
ling. Perhaps  the  most  significant  feature  is  that  super linear 
regions  occur  in  all  but  one  case.  One  possible  explanation  is  that 
these  regions  are  caused  by  space-charge-limited  (SCL)  current  flow 
arising  from  single  carrier  current  injection. 

Current  injection  phenomena  in  solids  have  been  reviewed  by 
Lampert  and  Mark  (Ref  10).  We  shall  outline  those  features  that  may 
be  relevant  to  this  work.  The  first  criterion  for  injection  to  occur 
is  the  existence  of  an  injecting  contact,  one  that  by  definition 
supplies  a reservior  of  carriers  freely  available  to  enter  the  insula- 
tor as  needed.  Lead  azide  has  a band  gap  of  approximately  3.5 eV 
(Ref  11,  12)  and  the  electron  affinity  has  been  found  to  be  about 
1.2eV  (Ref  13).  With  gold  contacts  on  lead  azide  single  crystals, 
it  is  expected  that  near  the  azide-electrode  interface,  the  lead 
azide  valence  band  may  be  bent  to  the  Fermi  level,  forming  an  inject- 
ing contact  for  holes.  This  expected  energy  band  scheme  is  schema- 
tically depicted  in  Figure  13a.  Because  of  its  relatively  large 
work  function,  gold  would  not  be  expected  to  form  an  injecting  con- 
tact for  electrons.  It  has  been  shown  (Ref  10)  that  in  the  presence 
of  thermally  generated  free  carriers  and  no  carrier  traps,  SCL  cur- 
rent flow  is  characterized  by  an  ohmic  I-V  relation  at  low  voltages, 
and  a quadratic  one  at  higher  voltages.  The  presence  of  traps  can 
lead  to  discontinuities  in  the  IV  characteristics. 
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The  problem  with  the  SCL  current  interpretation  is  that  recent 
internal  photoemission  studies  (Ref  14)  indicate  that  gold  forms  a 
blocking  (non-injecting)  contact  to  PbN^  under  the  preparation  condi- 
tions used  here.  This  is  depicted  in  Figure  13b.  Thus  the  super- 
linear  I-V  relations  cannot  be  explained  by  SCL  current  flow.  Either 
high-field  effects  (in  which,  e.g.,  the  mobility  is  a function  of 
field),  or  field  lowering  of  the  (blocking)  barrier  might  then  be 
responsible  (Ref  15). 

It  is  not  possible  at  this  time  to  identify  the  phenomena  which 
lead  to  initiation.  Some  comments  are,  however,  possible.  Initiation 
is  not  caused  by  bulk  joule  heating.  The  power  input  prior  to  initia- 
tion, IV,  is  typically  10" 5 watts.  Total  energy  input  up  to  initia- 
tion is  of  the  order  of  10" 3 calories  distributed  over  a volume  of 
10"2cm3  (.2cm  x .04cm2).  For  lead  azide,  the  density  (Ref  16)  is 
4.7gm/cm^  and  the  specific  heat  (Ref  1)  0. lOcal/gm-K.  Even  when  all 
energy  loss  processes  are  neglected,  this  leads  to  temperature  in- 
creases of  only  about  0.2K,  which  would  not  be  expected  to  cause 
initiation. 

Temperature  rises  could  be  much  higher  if  the  current  density 
were  concentrated  near  one  electrode,  for  example,  or  in  a small- 
diameter  filament  stretching  from  one  electrode  to  the  other.  The 
first  possibility  can  be  rejected  partly  because  the  relatively  con- 
stant threshold  electric  field  value  suggests  that  the  voltage  drops 
rather  uniformly  across  the  bulk  of  the  crystal.  Further,  the  rapid 
rate  of  thermal  diffusion,  in  this  essentially  one -dimensional  geome- 
try, makes  it  impossible  to  achieve  thermal  initiation  by  joule  heat- 
ing no  matter  how  thin  a laminar  region  the  voltage  drop  occurs  over 
on  our  time  scale.  Indeed,  a temperature  rise  of  only  0.1K  is 
expected  (see  Appendix).  For  the  case  of  filamentary  conduction, 
thermal  initiation  by  joule  heating  is  also  impossible,  with  a 
temperature  rise  of  about  IK  (see  Appendix).  Finally,  for  the  case 
of  all  the  energy  input  being  confined  to  a small  point- like  region 
in  the  crystal  (a  rather  unusual  physical  situation),  thermal  ini- 
tiation is  indeed  possible  on  a time  scale  of  approximately  10“8  secs 
(see  Appendix) . Such  short  time  scales  for  initiation  run  counter  to 
our  observations  and  are  not  consistent  with  our  observed  thresholds. 

Other  possible  phenomena  that  could  lead  to  initiation  would 
appear  to  be  essentially  non- thermal.  Dielectric  breakdown  at  these 
low  current  densities,  for  example,  involves  a manyfold  multiplica- 
tion of  current  carriers  (electrons  and/or  holes).  Such  multiplica- 
tion requires  quanta  of  energy  in  the  electron-volt  range,  which  are 
in  extremely  short  supply  at  normal  temperatures.  Any  temperature 
rises  associated  with  exothermic  reaction  would  themselves  have  been 
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induced  by  the  presence  of  electron  and/or  hole  densities  in  excess 
of  the  small  density  available  in  the  absence  of  the  electric  field. 

McLaren  and  Rogers  (Ref  17)  have  reported  that  AgNj  initiates 
under  applied  fields  of  ca  250  V/cm.  Bowden  and  McLaren  (Ref  2)  sug- 
gest that  this  is  caused  by  field  emission  from  the  cathode.  In  this 
model  electrons  enter  the  crystal  with  sufficient  energy  to  produce 
decomposition  (because  of  space-charge  fields),  followed  by  self- 
heating and  explosion.  While  such  a model  may  possibly  apply  to  our 
results  in  PbNg,  two  factors  should  be  noted.  Our  currents  are  four 
orders  of  magnitude  lower  than  those  observed  in  AgNj  even  for  much 
higher  fields.  And,  as  we  have  mentioned,  the  voltage  drop  in  our 
experiments  appears  to  be  relatively  uniform  across  our  crystals  at 
least  up  to  or  near  initiation. 

Whereas  the  present  report  finds  that  average  fields  of  3.5  x 
104  V/cm  lead  to  initiation  for  PbN^  crystals  using  Au-Au  contacts, 
PbN  crystals  and  pellets  with  capacitively  coupled  electrodes  can 
withstand  average  fields  of  1-2  x 105  V/cm  without  visible  damage 
(Ref  18) . The  role  of  the  contacts  is  thus  seen  to  be  important  in 
field  initiation,  and  a combined  experimental  and  theoretical  effort 
to  examine  that  role  will  be  reported  elsewhere  (Ref  18). 
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Table  1 


Electric  field  effects  on  lead  azide 


Sample 

Thickness 

L(cm) 

Current  at 
detonation 
(amp) 

Voltage  at 
detonation 

VDET 

3-72-30 

.019 

1.2 

X 

10-3 

800 

3-72-31 

.024 

6.9 

X 

10"  8 

1000 

3-72-32 

.022 

7.0 

X 

10-9 

800 

3-72-33 

.043 

1.6 

X 

h-» 

o 

1 

00 

1300 

3-72-34 

.025 

1.4 

X 

10"8 

1000 

3-72-35 

.076 

5.5 

X 

19"9 

2700 

3-72-36 

.073 

1.1 

X 

10"  8 

2400 

Threshold  field 
det/L  (volts/cm) 

4.20  x 104 
4.17  x 104 
3.60  x 104 
3.02  x 104 
4.00  x 104 
3.54  x 104 
3.26  x 104 
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Fig  1 Single  crystal  of  PbNg  with  evaporated  gold 
electrode.  (Dimension  ^0.3  x 0.4  x .02  cm) 


SCHEMATIC  OF  SAMPLE  MOUNTING 


CRYSTAL 


Fig  2 Crystal  mounted  for  measurements 


SAMPLE 


PICOAMMETER 


Fig  3 Schematic  of  measuring  circuit 


Fig  4 Sample  chamber 
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Fig  5 Typical  current  response  to  step 
Increases  in  voltages 
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Fig  6 Current  - Voltage  characteristic 
(.019  cm  thick) 
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Fig  7 Current  - Voltage  characteristic 
(.024  cm  thick) 
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CURRENT  (AMPS) 


Fig  8 Current  - Voltage  characteristic 
(.022  cm  thick) 
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CURRENT  (AMPS) 


Fig  10  Current  - Voltage  characteristic 
(.025  cm  thick) 
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CURRENT  (AMPS) 


Fig  11  Current  - Voltage  characteristic 
(.076  cm  thick) 
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Fig  12  Current  - Voltage  characteristic 
(.073  cm  thick) 
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Fig  13  Band  picture  of  (a)  an  ohmic  contact  for  holes 
and  (b)  a blocking  contact  for  holes;  Ep  is  the 
Fermi  level,  and  the  barrier  height 
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It  was  shown  in  the  text  that  for  our  experiments,  bulk  joule 
heating  of  the  entire  crystal  would  lead  to  insignificant  temperature 
rises.  We  examine  here  the  consequences  of  more  localized  joule  heat- 
ing, as  would  occur  when  either  the  current  or  voltage  drop  is  not 
uniformly  distributed.  Some  conclusions  are  possible  on  the  basis  of 
rather  simple  models.  For  purposes  of  comparison,  the  explosion 
temperature  (minimum  temperature  for  ignition  in  5 seconds)  for  lead 
azide  is  315°C  (Ref  1). 

Our  basic  assumptions  are:  Energy  input  occurs  primarily  through 

joule  heat,  i.e.,  self-heating  by  exothermic  chemical  decomposition  is 
not  important  at  temperatures  well  below  the  explosion  temperature. 
Energy  delocalization  is  effected  primarily  by  thermal  conduction. 

This  is  certainly  valid  for  a lower  limit  (upper  limit  on  temperature 
rise),  as  other  modes  of  energy  transport  (such  as  radiative  process- 
es) can  only  increase  the  rate  of  energy  delocalization.  The  material 
properties  of  lead  azide  needed  are  the  density  p and  the  specific 
heat  c (values  appear  in  the  text),  and  the  thermal  conductivity  k. 

A value  of  4.0  x 10’1*  cal/°C  cm-sec  is  assumed  for  k (see  Ref  1, 
although  there  is  ambiguity  in  the  units  quoted).  All  are  assumed 
to  be  independent  of  temperature.  The  value  for  k is  for  a pressed 
pellet.  One  would  expect  a higher  value  for  single  crystals,  so  that 
again  the  assumptions  are  consistent  with  a lower  limit  for  rate  of 
energy  delocalization.  The  quantities  p and  c will  appear  in  the 

thermal  diffusivity  k,  defined  as  k = — . 


One- Dimensional  Geometry 

We  first  examine  the  temperature  rise  with  time  for  the  case 
where  the  voltage  drop  occurs  across  an  infinitely  thin  plane  paral- 
lel to  the  electrodes.  This  can  happen  either  within  the  lead  azide 
crystal  or  near  one  of  the  electrodes.  The  physical  system  is  then 
reasonably  approximated,  at  least  for  an  upper  limit  of  temperature 
rise  t,  by  the  adjoining  semi-infinite  regions  depicted  in  Figure  A-l 
with  thermal  diffusivity  k and  thermal  conductivity  k defined  in  each 
region.  For  time  t>0,  heat  is  supplied  at  the  constant  rate  FQ  per 
unit  time  per  unit  area  in  the  plane  x=0. 

The  thermal  flux  f is  then  defined  as 


f = -k  -21  (A-l) 

at 
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and  satisfies  the  differential  equation 


< a2f  = 

2 3x2  9t 


where  Kj 
The 


is  used  for  x<0  and  f°r  ^>0. 

solution  for  t,  by  well-known  methods  (Ref  19),  is 


T = 


2Fi  “i  l) 


1/2 


ierfc 


2(Kl  t) 
2 


1/2' 


(A-2) 


The  upper  subscript  is  used  for  x<0  and  the  lower  for  x>0.  The  func 
tion  ierfc  (A)  is  defined  as 


ierfc  (A)  = — e~*  - A erfc  (A) , (A-3) 

/iT 

where 

QO 

erfc  (A)  = / e dy.  (A-4) 

^ A 

The  latter  function,  commonly  occurring  in  heat  conduction  problems, 
is  tabulated  in  a number  of  placed  (Ref  19) . The  quantities  Fj  and 
?2  are  determined  by  the  conditions 


1/2 

1 K 


k 


1/2 


2 k2 


and  + ^2  = 


(A— 5) 


The  highest  temperature  at  any  given  time  occurs  at  x=o,  where 
the  temperature  tq  is 


= 1.128  F2  (*2t) 


1/2 


(A- 6) 
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If  the  sudden  voltage  drop  occurs  well  within  the  crystal,  clearly 
Fj  = F0  = FQ/2.  Under  the  conditions  of  our  experiment,  FQ  6 x 

10~5cal/sec-cm2  when  the  voltage  is  on.  Plugging  in  values  for  the 
thermal  properties  k and  k,  one  finds  (for  t in  seconds)  t0  * 2 x 

10-3  t1/2K.  Thus  for  t = 103  secs,  the  maximum  temperature  rise 
resulting  from  joule  heating  for  this  case  of  rapid  voltage  drop  is 
less  than  0.1K.  Should  the  sudden  voltage  drop  occur  near  a gold 
contact,  the  corresponding  temperature  rise  calculated  using  Equa- 
tions (A-5)  and  (A-6)  is  two  orders  of  magnitude  lower  because  of 
the  high  thermal  conductivity  of  gold.  We  conclude  that  even  for 
abrupt  voltage  drops,  our  experimental  power  input  is  insufficient 
to  cause  thermal  initiation  of  lead  azide  by  simple  joule  heating 
temperature  rises. 


Two-Dimensional  Geometry 

We  next  examine  the  case  of  all  the  current  flowing  in  a small 
filament  joining,  and  perpendicular  to,  the  electrodes.  A reasonable 
idealization  of  this  case  is  an  infinitely  thin  constant  line  source 
of  heat  which  begins  emitting  at  t=0,  embedded  in  an  infinite  lead 
azide  crystal. 

Radial  distance  from  the  line  source  is  denoted  by  r.  For  this 
case,  we  search  first  for  a solution  to  the  heat  conduction  differ- 
ential equation 


V2t  = I !i. 
k 3t 


(A- 7) 


subject  to  release  at  t=0  of  an  instantaneous  uniform  heat  pulse  of 
magnitude  Q pc  (p  is  density,  and  c specific  heat)  per  unit  length 
of  the  line.  The  solution  for  all  times  t>0  can  be  shown  to  be 


x 


— e“r  /4<t 
4lTKt 


(A- 8) 


As  t-+0,  t also  goes  to  zero  for  all  r except  r=0.  There  it  be- 
comes infinite  in  such  a way  that  its  integral  over  all  space  is 
finite.  Physically  this  occurs  because  a finite  amount  of  energy 
was  assumed  supplied  to  an  infinitesimally  small  volume,  that  of  the 
ideal  line.  In  practice,  heat  cannot  be  released  from  an  ideal 
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line  source.  Indeed,  temperature  is  defined  only  in  a statistical 
sense,  and  it  therefore  has  no  physical  meaning  when  applied  to  re- 
gions of  sub-atomic  dimensions,  i.e.,  for  r<10"8cm. 


If  now  heat  is  released  by  the  source  continuously  at  the  con- 
stant rate  pcq  per  unit  time  per  unit  length,  the  temperature  is 
determined  by  integrating  Equation  (A-8)  over  time,  i.e.. 


T 


T 

a3"  1 
4tt<  0 


dt 

i 

t 


(A-9) 


Integration  yields  the  function 


T 


= . _SL 

4ttk 


Ei 


(A- 10) 


where 


-Ei  (-x) 


I 

X 


u 


For  small  values  of  x, 


Ei  (-x) 


Y + ln(x) 


X + + 0(x3)  , 


(A-ll) 


where  y = 0.5772.  For  our  present  purposes,  it  will  be  the  case  that 
Ei(-x)  is  sufficiently  approximated  by  the  single  term  ln(x). 

Under  the  conditions  of  our  experiment  pcq  * 10-1*  cal/sec-cm. 

For  r = 10'8  cm  and  t = 103  sec,  one  finds  from  Equation  (A- 10)  that 
t<1K.  Temperature  rises  of  the  order  100K  can  be  realized  for 
r 88  io" 1000  cm,  but  we  have  shown  that  the  concept  of  temperature  is 
not  physically  meaningful  for  these  dimensions.  We  conclude  that 
even  for  the  case  of  all  the  current  flowing  in  a narrow  filament 
through  the  crystal,  our  experimental  power  input  is  insufficient  to 
cause  thermal  initiation  of  lead  azide  by  simple  joule  heating  tem- 
perature rises. 
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Three-Dimensional  Geometry 

We  finally  consider  the  case  of  both  filamentary  conduction  and 
an  abrupt  voltage  drop,  so  that  all  the  energy  input  is  delivered  in 
a small  volume  of  the  crystal.  This  situation  can  be  approximated 
by  a continuous  point  source  in  an  infinite  medium. 

The  solution  to  the  heat  conduction  differential  Equation  (A- 7) 
for  an  instantaneous  heat  point  source  of  magnitude  pcP,  located  at 
r=0  and  delivered  at  t=0,  can  be  shown  to  be 


T 


P 

(4ir<t)  3^2 


-r2 
4<t  ’ 


(A- 12) 


for  times  t>0.  The  solution  for  a continuous  steady  heat  point  source 
of  magnitude  p per  unit  time  is  then  obtained  by  integrating  Equa- 
tion (A- 12)  over  time,  and  is 


t = P erfc  I 1/2. 

emcr  (4<t) 


(A- 13) 


The  function  erfc  (A)  has  been  defined  in  Equation  (A-4) . 

p 

For  t-K»  or  for  sufficiently  small  r,  this  reduces  to  t = ^Tiricr' 
Under  the  conditions  of  our  experiments  p as  2 x 10-6  cal/sec.  It 
follows  that  temperature  rises  in  excess  of  100K  are  realized  within 
a radius  of  r « 4 x 10~6cm.  We  conclude  that  initiation  by  simple 
joule  heating  temperature  rises  is  possible  if  all  the  power  input 
is  confined  to  a small  point-like  region  of  the  crystal.  Under  these 
conditions  initiation  should  in  fact  occur  on  a time  scale  of  approxi- 
mately 10-8  sec.  The  fact  that  our  experiments  show  initiation  occur- 
ring after  many  seconds  tends  to  discredit  the  applicability  of  this 
type  of  mechanism.  Further,  this  kind  of  mechanism  would  predict 
initiation  for  far  lower  fields  than  those  observed  to  be  necessary. 
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Fig  A1  One-dimensional  model  for  voltage  drop  across 
a thin  plane  at  x=0.  The  adjoining  media  are 
characterized  by  thermal  diffusivities  and 
thermal  conductivities  , k1  and  <2>  k,,. 
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